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We investigate the phase structure of lattice QCD for general number of flavors Np. 
Based on numerical results combined with the result of the perturbation theory we 
propose the following picture; When Np > 17, there is only a trivial fixed point 
and therefore the theory in the continuum limit is trivial. On the other hand, when 
16 > Np > 7, there is a non-trivial fixed point and therefore the theory is non- 
trivial with anomalous dimensions, however, without quark confinement. Theories 
which satisfy both quark confinement and spontaneous chiral symmetry breaking 
in the continuum limit exist only for A'^^' < 6. 



1 Introduction 

The fundamental properties of QCD are quark confinement, asymptotic free- 
dom and spontaneous breakdown of chiral symmetry. It is well-known that if 
the number of flavors exceeds 16^, the asymptotic freedom is lost. Then, the 
question which naturally arises is whether there is constraint on the number 
of flavors for quark conflnement and/or the spontaneous breakdown of chiral 
symmetry. 

In our previous workld, it was shown that in the strong coupling limit, when 
the number of flavors Np is greater than or equal to 7, quarks are deconfined 
and chiral symmetry is restored at zero temperature, if the quark mass is 
lighter than a critical value. In this work we investigate the problem of what is 
the condition on the number of flavors for quark confinement and spontaneous 
chiral symmetry breaking in the continuum limit □. 

As asymptotic freedom is the nature at short distance, one can apply the 
perturbation theory to investigate the critical number for it. However, because 
quark confinement and spontaneous chiral symmetry breaking are due to non- 
perturbative effects, one has to apply a non-perturbative method throughout 
the investigation of the critical numbers for them. Therefore we employ lattice 
QCD for the study in this work, since lattice QCD is such a theory constructed 
non-perturbatively. 

In order to investigate the condition on the number of flavors Np in the 
continuum limit, we first clarify the phase structure at zero temperature for 
general number of fiavors, in particular, for Np > 7. When the phase diagram 
becomes clear, we are able to see what kind of continuum limit exists and 
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eventually answer the question given above. 

In Sec. 1^ the formulation of lattice QCD is introduced. The notations 
which will be used later, and the Wilson quark action and the standard one- 
plaquette gauge action which we employ in this work are given. In Sec. I 
quark mass is defined and chiral property in the Wilson formalism of fermions 
is discussed. Before going into details a brief summary is given Sec. ^. In 
Secs.^ - |l^, details of numerical simulations and results are given. Conclusions 
are given m Sec. 0. 

2 Lattice QCD 

Lattice QCD is defined on a hypercubic lattice in 4-dimensional euclidean space 
with lattice spacing a. A site is denoted by a vector n = (ni, ^2,^3, 724), where 
rii's are integers. A link with end points at the sites n and n + fi is specified 
by a pair {n,fj,), where fi denotes a unit vector in the /i direction. 

The gauge variable J7„,ai which is an element of SU(3) gauge group is 
defined on the link (n, /i). The action for gluons which we adopt is given by 



where g is the gauge coupling constant. This action is called the standard one- 
plaquette gauge action. We usually use, instead of the bare gauge coupling 
constant g, (3 defined by 



The quark variable is a Grassman number defined at a site. It is well- 
known that a naive discretization of the Dirac action 



with Too being the bare fermion mass, leads to 16 poles instead of one pole. 
This problem is called the species doubling. To avoid this problem K. Wilson 
proposed adding to the naive discretized action a dimension 5 operator called 
the Wilson term 




(1) 




(3) 




(4) 
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Rescaling the fermion field by q = \/2Kip and making the action gauge invari- 
ant, we obtain the Wilson fermion action 

n 

where 

''=(moa'+4r)' 

which is called the hopping parameter. 

The full action S is given by the sum of the gauge part Sgauge and the 
fermion part SwiUon, 

S = Sgauge + S\Yilson- (7) 

The expectation value of an operator 0{U, ip, 4') is given by 

<0>=^ In dUn,^ n di;nd^nO{U, i>)eMS), (8) 

where Z is the partition function 

Z= In d[/„,^ n dV'„dV'„exp(5) , (9) 

n,iJ, n 

with dUn^^ being the Haar measure of SU(3). 

In the case of degenerate Np flavors, lattice QCD contains two parameters: 
the gauge coupling constant (3 = 6/5^ and the hopping parameter K. In the 
non-degenerate case, the number of the hopping parameters is Np- 

The lattice spacing o is a function of the bare coupling constant which is 
defined through the RG beta function. The continuum limit is the limit where 
the lattice spacing a goes to zero in such a way that physical quantities are 
kept constant. In an asymptotically free case the limit corresponds to the limit 
where the coupling constant g vanishes. 

At finite temperatures the linear extension in the time direction iVj is much 
smaller that those in the spatial directions {Nx,Ny,Nz). The temperature T 
is given by 

l/Nta. 

For gluons the periodic boundary condition is imposed, while for quarks an 
anti-periodic boundary condition is imposed in the time direction. 
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3 Quark mass and chiral symmetry in the Wilson quark action 

In the formalism of Wilson of quarks on the lattice, the flavor symmetry as well 
as C, P and T symmetries are exactly satisfied on a lattice with a finite lattice 
spacing. However, chiral symmetry is explicitly broken by the Wilson term 
even for the vanishing bare quark mass mo — 0. The lack of chiral symmetry 
causes much conceptual and technical difficulties in numerical simulations and 
physics interpretation of data. (See for more details Ref.13 and references cited 
there.) 

The chiral property of the Wilson fermion action was_first systematically 
investigated through Ward identities by Bochicchio et al. a. We also indepen- 
dently proposedu to define the current quark mass by 

2mq < I P I TT >= -TOtt < I I TT > 

where P is the pseudoscalar density and A4 the fourth component of the local 
axial vector current. We use this definition of the current quark mass in this 
work. In general we need a multicative normalization factor for the axial 
current, which we have absorbed into the definition of the quark mass, because 
this definition is sufficient for later use. The value of the quark mass does not 
depend on whether the system is in the high or the low temperature phase 
when the gauge coupling constant is not so large: /3 > 5.5 for Np — 2. 

With this definition of quark mass, the pion mass vanishes in the chiral 
limit where the quark mass vanishes at zero temperature, for Np < 6. However, 
in the deconfining phase at finite temperatures the pion mass does not vanish 
in the chiral limit. It is almost equal to twice the lowest Matsubara frequency 
n/N^. This implies that the pion state is a free two-quark state. Quarks are 
not confined. The pion mass is nearly equal to the scalar meson mass, and the 
rho meson mass to the axial vector meson mass. The chiral symmetry is also 
manifest within corrections due to finite lattice spacing. 

Thus hadron masses in the chiral limit where the current quark mass van- 
ishes clarifies the phase of the system at finite temperatures for Np < 6. We 
use the same strategy to discriminate the phase of QCD with various number 
of flavors Np: the pion mass at zero quark mass determines whether chiral 
symmetry is spontaneously broken or not, and also whether quarks are con- 
fined or not. It should be noted that in QCD with Np fiavors there is no known 
order parameter for quark confinement. When the pion mass is about twice 
the lowest Matsubara frequency we conclude that quarks are not confined. 
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Figure 1: Our conjecture for the beta fuction 

4 Brief Summary 

Before going into details let me give a brief summary, because the following 
sections will be too technical for non-experts of lattice QCD. 

In Fig. |l] is presented our conjecture in terms of the beta function of renor- 
malization group. Here g is the bare gauge coupling constant. When Np is 
equal or smaller than 6, the beta function is negative for all coupling constants. 
On the other hand, when Np is equal or larger than 17, it is positive. In the 
case when Np is between 7 and 16, it changes sign from negative to positive 
with increasing g. Therefore there is a non-trivial IR fixed pointy The existence 
of such a fixed point was first pointed out by Banks and Zaks.O However, our 
picture is diffcreiit those proposed by them. Fig. ^ is the conjecture made by 
Banks and ZaksQ. They assumed that the beta function becomes negative for 
large enough coupling constant for all Np. This point is different from ours. 
Their critical value is 9, which is also different from our value 7. 

It is sometimes argued in the literature that quarks are confined in the 
strong coupling constant. Although there is a proof for quark confinement in 
the strong coupling constant in the pure SU{3) gauge theory, there are no 
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Figure 2: Conjecture for the beta fuction by Banks and Zans 

such proofs in QCD with dynamical quarks. Our previous result is a counter 
example for such an argument. 

Fig. ^ shows the phase structure in the /3 - K plane for various number 
of flavors. When Np < 6, there is the chiral limit Kc{/3), where the current 
quark niq vanishes, in the confining phase. The value of Kc{P) at /3 = cxd is 1/8, 
which corresponds to the vanishing bare quark mass toq = 0. As (3 decreases, 
the Kc increases, up to 1/4 at /3 = 0. If the action would be chiral symmetric, 
the chiral line should be a constant, 1/8. The line K = corresponds to an 
infinite quark mass. Quarks are confined for any value of the current quark 
mass for all values of f3. 

When Np > 7, there is no chiral limit in the confining phase. There is a 
line of a first order phase transition from the confining phase to a deconfining 
phase at a finite current quark mass for all values of (3. It is crucial that the 
point X = 1/8 at /? = oo does belong to a deconfining phase. The chiral line 
which starts from this point hits the boundary between the confining phase 
and the deconfining phase at finite /3 as shown in Fig. ^ for the case where 
Np is not too large. This situation makes the interpretation of numerical data 
difficult; the data around (3 = 4.5 for Np = 7, 12 and 18 look very complicated. 
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Figure 3: Phase structure for various number of flavors Np and the RG flow. 



In order to understand the structure of the deconfining phase more clearly, 
we increase Np up to 300, because the region of the deconfining phase be- 
comes larger with larger Np, and consequently the structure becomes simple, 
as shown in Fig. |^. When Np > 240, the chiral line where the quark mass van- 
ishes goes straight from (3 — oo to /? = 0, without hitting the phase transition 
boundary. We have made a MCRG study to investigate the flow of RG along 
the massless line. The direction at /? = c» is known from the perturbation 
theory, from left to right in the /3 - K plane. Our results of a MCRG study 
which will be described later implies that the direction is identical all through 
the /3 region. Thus there is only a trivial IR fixed point at /3 = oo. There are 
no other fixed points where we are able to take a continuum limit of the theory 
with quarks. Thus the theory is a trivial free theory. 

As Np decreases, the massless quark line hits the boundary of a first order 
phase transition. Except this point other general features are the same as 
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Figure 4: The phase structure for Np < 6 at zero temperature and at finite temperatures. 



those for Np > 240. Massless quarks exist only in the deconfining phase. A 
continuum hmit can be taken only at a fixed point on this line. Although we 
have not investigated the flow of RG, we naturally assume that the direction 
for Np < 17 is the same as that for Np > 240. Thus the theory is also trivial 
for Np< 17. 

When 7 < Np < 16, the phase structure seems to be not too much different 
from those for Np > 17. However, the flow of RG at /3 = oc' is opposite 
to that for Np > 17, because the theory is asymptotically free. Therefore a 
continuum limit is governed by an IR fixed point which should exist somewhere 
at finite value of /3. Unfortunately, as explained above, the massless line hits the 
boundary in the two parameter space. Therefore it is hard to see numerically 
where the IR fixed point exists. When one makes a RG transformation the 
trajectory moves into, in general, a multi-parameter space. Thus within the 
two parameter space we are investigating here it is difficult to pin down the 
position of the IR fixed point. However, even without knowing the position 
we are able to discuss the nature of the continuum limit. On the massless line 
around (3 = 4.5, the mass of the pion is roughly twice the lowest Matsubara 
frequency. This implies that the quark is almost free. That is, an anomalous 
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dimension which is governed by the IR fixed point is smaU. The fact that the 
anomalous dimension is smaU suggests that the IR fixed point exists at finite (3: 
If the IR fixed point would exist at = or 5 = cx), the anomalous dimension 
would be large or infinite. Anyway, the fact that the mass of the pion is roughly 
twice the lowest Matsubara frequency implies that quarks are not confined in 
this phase and the chiral symmetry is not spontaneously broken. Thus this 
phase is a phase governed by a non-trivial IR fixed point. 

A salient fact for the above in this section is that in the strong coupling 
limit, when the number of flavors Np is greater than or equal to 7, quarks are 
deconfined and chiral symmetry is not spontaneously broken at zero tempera- 
ture, if the quark mass is lighter than a critical value, which was shown in our 
previous workl^. 

In the following sections details of numerical simulations will be given. 
One has to perform numerical simulations on a lattice with a finite Nt. This 
implies that when (3 becomes large one encounters a finite temperature phase 
transition. A schematic diagram is shown in Fig. ^ for the case Np < 6. 
In the case Np > 7, one has to be careful about the existence of both zero 
temperature transitions and finite temperature transitions. 

5 Simulation parameters 

The lattice sizes are 8^ x 10 x Nt (Nt = 4, 6 or 8), 16^ x 24 x iVf {Nt ^ 16) 
and 18^ x 24x Nt {Nt = 18). We use an anti-periodic boundary condition for 
quarks in the t direction and periodic boundary conditions otherwise. When 
the hadron spectrum is calculated, the lattice is duplicated in the direction 
of lattice size 10 for Nt < 8, which we call the z direction. We call the pion 
screening mass simply the pion mass..and similarly for the quark mass. 

We use the hybrid R algorithm Ll for the gauge configuration generation. 
As Np increases we have to decrease Ar, such as Ar =0.0025 for Np = 240, to 
reduce O(At^) errors. We have checked that the errors are sufficiently small 
selecting typical cases. 

6 Kd versus Np 

We first study the case of /3 = as in the previous worki. Fig. |^a shows the 
results of the deconfining transition points Kd at /3 = obtained on A^f = 4 
lattices for various numbers of Np (see also Fig. ||d). 

The data of plaquette for Np — 240 at /3 = indicates that the decon- 
fining transition is of first order and that the location of the transition point 
is independent of Nt for large Nt {l/Kd ~ 8.1 for A^t = 4 and l/Kd ~ 7.75 
for Nt = 8 and 16). Therefore we conclude that the transition is bulk as is 
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1/K 



Figure 5: a) The transition point 1/Kd at /3 = versus Np at Nt = 4 and Nt > 8. For 
clarity, data at A^t = 8 for Np =1-800 is slightly shifted in the figure. Shaded region was 
investigated in our previous study u. b) Plaquette versus 1/K at /3 = for various number 
of Np at Aft = 4. 



confirmed in our previous work for Np = 7\^. In Fig. ^ the transition points 
at A^t > 8 for Np = 7, 240 and 300 are also included. These values are roughly 
those for the bulk transition points at zero temperature. 

As Np increases up to 240, 1/Kd approaches toward or exceeds the value 
8, which is the value of 1/K for a massless free quark. Because of this, for 
Np > 240, we are able to see a wide region of the deconfining phase in the 
whole range of /3. Therefore we intensively investigate the case Np — 240, and 
then decrease Np. 
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7 Nf = 240 

Fig. I shows the results of ml and 2mg for Np = 240 at /3 = 0, 2.0, 4.5, 6.0, 
and 100 on the Nt = A lattice. A very striking fact is that the shape of and 
2mq as a function of 1/i^ only slightly changes for 1 /K < 8 when the value of 
(3 decreases from oo down 0. (Note that the values at /3 = oo are those for free 
quarks.) Only the position of the local minimum of to^ at 1/K ~ 8, which 
corresponds to the vanishing point of m^, slightly shifts toward smaller 1/K. 
We obtain similar results also for Nt — 8. 

From these data we propose the phase diagram in Fig. I for Nf = 240. 
The dark shaded line is the phase boundary between the confining phase and 
the deconfining phase at zero temperature. When A^t = 4 or 8, the boundary 
line bends down at finite /3 as shown in Fig. |7| due to the finite temperature 
phase transition of the confining phase. 

The dashed line corresponds to the niq — line. This line also corresponds 
to the minimum point of mj. We have also calculated the quark propagator 
in the Landau gauge and checked that chiral symmetry of the propagator, 
75G(z)75 = — G(z), is actually satisfied on the iriq = line. 

The results for the case Np = 300 are essentially the same as those for 
Np ~ 240 except for very small shifts of the transition point and the minimum 
point of m^. 
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Figure 7: Phase diagram for Np = 240. 



8 Renormalization group flow 

The niq — point at /3 = cx) is the trivial IR fixed point for Np > 17. The 
phase diagram shown in Fig. ^ suggests that there are no other fixed points on 
the TTiq = hne at finite /3. 

In order to confirm this, we investigate the direction of Renormalization 
Group (RG) flow along the rriq — line for Np — 240, using a Monte Carlo 
Renormalization Group (MCRG) method. 

One problem here is that it is practically impossible to make simulations 
in the massless limit at zero temperature due to the existence of zero modes in 
the quark matrix. Therefore, we impose an anti-periodic boundary condition 
in the t direction. 

We make a block transformation for a change of scale factor 2, and estimate 
the quantity A/3 for the change of a a' = 2a. We generate configurations on 
an 8^ lattice on the rriq = points at /3 = and 6.0 and make twice blockings. 
We also generate configurations on a 4"^ lattice and make once a blocking. Then 
we calculate A/3 by matching the value of the plaquette at each step.^ 

From the matching, we obtain A/3 ~ 6.5 for /3 = and 10.5 for (3 — 6.0. 
The value obtained from the perturbation theory is A/3 ~ 8.8 at (3 — 6.0 for 
Np = 240. The signs are the same and the magnitudes are comparable. This 
suggests that the directions of RG flow on the iJiq = line at /3 = and 6.0 
are the same as that at /3 = cx) for Np — 240. This further suggests that there 
are no fixed points at finite /3. All of the above imply that the theory is trivial 

" It is known for the pure SU{3) gauge theory, in particular in the deconfining phase, that 
one has to make a more careful analysis using several types of Wilson loop to extract a 
precise value of A/3. We reserve elaboration of this point and a fine tuning oi 1/K at each 
/3 for future works. 
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Figure 8: Phase diagram for Np = 18. 



in the case of Np = 240. 
9 240 > A^i^ > 17 

Now we decrease the value of Np from 240. When (3 = 6.0, the shapes of 
are almost identical to each other, except for a slight shift toward smaller 1/K 
as Np decreases. 

When /3 = 0, the boundary of the first order phase transition between 
the deconfining phase and the confining phase moves toward smaller 1 / K and 
therefore the range of the deconfining phase decreases: 1/ Kd ~ 8.5, 8, 7.6, 7.2, 
6.1 for Np = 300, 240, 180, 120 and 60, respectively, as shown in Fig. |. For 
Np = 18, the value of 1/Kd decreases down to 5.0. 

Due to the fact that the confining phase invades the deconfining phase, 
the massless line in the deconfining phase hits the boundary at finite [3 when 
Np becomes small. For example, in the case of Np = 18, as shown in Fig. |[ 
it hits at /3 = 4.0 — 4.5. 

Although the area of the deconfining phase decreases with decreasing Np, 
the shape and the position of in the part of deconfining phase only slightly 
change from Np ~ 300 to 18. The values of ruq as functions of 1/K also show 
only slight changes in the deconfining phase. These facts, combined with the 
perturbative result that for Np > 17, /3 = oo is the IR fixed point, suggest that 
the RG flow along the massless quark in the deconfining phase is the same as 
that for Np = 240. In this case, the theory is trivial for Np > 17.[| 

Although we have investigated only down to the Np = 18 case, we do not expect any 
qualitative differences between the cases of Np = 17 and 18. 
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Figure 9: Phase diagram for Np = 12. 
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10 16 > > 7 

When Np < 16, the theory is asymptotic free. On the other hand, quatk 
confinement is lost for Np > 7 even in the strong couphng hmit /? = Id. 
Therefore the question is what happens for the cases 16 > Np > 7 in the 
continuum hmit. 

We have intensively simulated the cases Np ~ 12 and 7. The phase di- 
agram for Np = 12 is shown in Fig. |[ Although the gross feature of the 
phase diagram seems to be not too much different from the case of Np = 18 
which is shown in Fig. |8[ the flow of RG at f3 — oo is opposite to that for 
Np > 17, because the theory is asymptotically free. Therefore a continuum 
limit is governed by an IR fixed point which should exist somewhere at fi- 
nite value of (3. Unfortunately, as explained above, the massless line hits the 
boundary in the two parameter space. Therefore it is hard to see numerically 
where the IR fixed point exists. When one makes a RG transformation the 
trajectory moves into, in general, a multi-parameter space. Thus within the 
two parameter space we are investigating here it is difficult to pin down the 
position of the IR fixed point. Even without knowing the position we are able 
to discuss the nature of the continuum limit. On the massless line around 
f3 = 4.5, the mass of the pion is roughly twice the lowest Matsubara frequency. 
This implies that the quark is almost free. That is an anomalous dimension 
which is governed by the IR fixed point is small. The fact that the anomalous 
dimension is small suggests that the IR fixed point exists at finite /?: If the IR 
fixed point would exist at /? = or 5 = cx), the anomalous dimension would 
be large or infinite. Anyway, the fact that the mass of the pion is roughly 
twice the lowest Matsubara frequency implies that quarks are not confined in 
this phase and the chiral symmetry is not spontaneously broken. Thus this 
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phase is a phase governed by a non-trivial IR fixed point. The theory in the 

continuum limit is a non-trivial theory with anomalous dimensions, however, 
without confinement. We reserve for future works a detailed RG study to find 
a non-trivial fixed point in a wider parameter space. 

11 Conclusions 

Based on numerical results combined with the result of the perturbation theory 
we propose the following picture: There are three categories depending on 
the number of flavors Np: a free theory for Np > 17, a non-trivial theory 
governed by a non-trivial IR fixed point for 16 > Np > 7 and confinement and 
spontaneous chiral symmetry breaking for Np < 6. 
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